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A first-principles theory of resonant magnetic scattering of x rays is presented. The scattering 
amplitudes are calculated using a standard time-dependent perturbation theory to second order 
in the electron-photon interaction vertex. In order to calculate the cross section reliably an ac- 
curate description of the electronic states in the material under investigation is required and this 
is provided by the density functional theory (DFT) employing the Local Spin Density Approx- 
imation combined with the self-interaction corrections (SIC-LSD). The magnetic x-ray resonant 
scattering (MXRS) theory has been implemented in the framework of the relativistic spin-polarized 
LMTO-ASA band structure calculation method. The theory is illustrated with an application to 
ferromagnetic praseodymium. It is shown that the theory quantitatively reproduces the dependence 
on the spin and orbital magnetic moments originally predicted qualitatively (Blume, J. Appl. Phys, 
57, 3615 (1985)) and yields results that can be compared directly with experiment. 

PACS numbers: PACS numbers: 78.70.Ck, 75.25. +z, 71. 15.-m,71.20.Eh 



I. INTRODUCTION 



Magnetic X-Ray Scattering (MXRS) is a well devel- 
oped technique for probing the magnetic and electronic 
structures of materials. The foundations of the theory of 
MXRS were laid down by Blume. Later on Blume and 
Gibbs developed the theory further to show that the 
orbital and spin contributions to the magnetic moment 
can be measured separately using MXRS with a judicious 
choice of experimental geometry and polarization of the 
x rays. Hannon et al. |3( presented a nonrelativistic 
theory of x-ray resonance exchange scattering and wrote 
down explicit expressions for the electric dipole (El) and 
quadrupole (E2) contributions. This work is based on 
an atomic model of magnetism and has been applied 
successfully to a variety of materials including UAs and 
Gd by Fasolino et al. Q Rennert [j| produced a semi- 
relativistic theory of MXRS written in terms of Green's 
functions, but no such calculations have been performed. 
More recently, theory based on an atomic model of the 
electronic structure of materials has been written down 
by Lovesey |6j and co-workers and applied successfully to 
a variety of materials. Takahashi et al. have reported a 
theory which includes the band structure in the calcula- 
tion of anomalous x-ray scattering. A first-principles 
theory of MXRS based on a time-dependent second order 
perturbation theory and densit y fu nctional theory 0, Q 
was produced by Arola et al. [ToL ITU and applied suc- 
cessfully to several transition metal materials. This 
theory is restricted in its range of application because of 
the limitations imposed by the local density approxima- 
tion to DFT which means that the theory can only be 
applied to simple and transition metal materials. This 
is particularly unfortunate because it is in the rare earth 
and actinide materials that the most exotic magnetism 



in the periodic table occurs. 

In recent years advances in electronic structure cal- 
culations beyond the local density approximation have 
broadened the range of materials for which numerically 
accurate electronic structure calculations can be per- 
formed. In particular the LDA+U method [13| and the 
self-interaction corrected local s pin de nsity ap proxima- 
tion to density functional theory [lj, E3 have met 

with considerable success in describing materials with lo- 
calized electrons. The latter method reduces the degen- 
eracy of the / states at the Fermi level and hence also 
circumvents all the convergence problems associated with 
the LSD approximation to DFT in electronic structure 
calculations for rare earth materials. Notably, the SIC- 
LSD has provided a very good description of the rare 
earth metal and rare earth chalcogenide crystal struc- 
tures. A relativistic version of the SIC formalism 
has been derived 0] that has been shown to yield an 
excellent description of the electronic structure of rare 
earth materials in the few cases to which it has been ap- 
plied. This method is reviewed by Temmerman et at. 

m 

The fact that electromagnetic radiation can be scat- 
tered from the magnetic moments of spin- 1/2 particles 
was first shown by Low, and by Gell-Mann and Gold- 
berger half a century ago. [2(j Later on it was Platz- 
man and Tzoar [2l| who first proposed the use of x-ray 
scattering techniques to study the magnetization density 
of solids. At that time progress in studying magnetic 
structures using x rays was severely hampered because 
the cross section for magnetic scattering is smaller than 
the cross section for charge scattering [j| by a factor of 
(hw/mc 2 ) 2 . It was Gibbs et al. |22| who first observed 
a large resonant enhancement of the cross section when 
the energy of the x ray is tuned through an absorption 



2 



edge. Since that time technological advances have pro- 
duced high resolution, high intensity synchrotron radia- 
tion sources that have transformed magnetic x-ray res- 
onant scattering into a practical tool for investigating 
magnetic and electronic structures of materials. Nowa- 
days the world's leading synchrotron facilities have beam- 
lines dedicated to this technique |23j and applications 
of resonant x-ray scattering are burgeoning. Reviews of 
the experimental state-of-the-art MXRS techniques have 
been written by Stirling [24j and Cooper. 

Other approaches to interpreting MXRS spectra exist, 
particularly the successful methods based on group the- 
ory and angular momentum algebra that result in sum 
rules as described by Borgatti [2|j and by Carra and 
Luo. [2^ The present work should not be regarded as 
a rival theory to these, but rather as an attempt to ex- 
tend the range of density functional methods to describe 
magnetic scattering of x rays in the same way as is done 
for photoemission and other spectroscopies. [22| As a 
DFT-based theory our work is, of course, based on very 
different approximations to this earlier work, making di- 
rect comparison between the two theories problematic. 

We have recently implemented a first-principles theory 
of MXRS that is based on a standard time-dependent 
perturbation theory where the scattering amplitudes are 
calculated to second order in the electron-photon inter- 
action vertex. To describe MXRS from a given material 
it is necessary to have an accurate description of the elec- 
tronic structure of the material in question. This is pro- 
vided by using the SIC within the LSD approximation to 
the density functional theory which is implemented using 
the relativistic spin-polarized LMTO-ASA band struc- 
ture calculation method. |3Tj The theory of MXRS is 
equivalent to that of Arola el al., |10| but has been rewrit- 
ten in a form that is appropriate for implementation in 
connection with the LMTO-ASA method where there is 
substantial experience of SIC methods. The major step 
forward reported in this paper is the integration of the 
SIC into the MXRS theory which enables us to describe 
rare earth and actinide materials on an equal footing with 
transition and simple materials. 

In this paper, we give a detailed description of the 
MXRS theory and illustrate it in a calculation for 
praseodymium. The results are analysed and discussed. 
Finally we show that the present work is consistent with 
the earlier theory and demonstrate how the MXRS cross 
section reflects the properties of these materials. 



proximation to the total energy functional, which leads to 
a greatly improved description of static Coulomb correla- 
tion effects over the LSD approximation. This has been 
demonstrated in studies of the Hubbard model, [3^, 01 
in applications to 3d monoxide s ITH (l7| and cuprates, 

jj. l3fij /-electron systems, 

381 metal-insulator transitions 




orbital ordering, 
and solid hydrogen. 



For many applications it is necessary to account for all 
relativistic effects including spin-orbit coupling in an elec- 
tronic structure calculation. Relativistic effects become 
progressively more important as we proceed to heavier 
elements. They are also extremely important when we 
are considering properties dependent on orbital moments 
and their coupling to electron spins. 

The relativistic total energy functional in the local spin 
density approximation is 



E LSD [n{r)} = E kin [n(r)] + U[n(r)]+ I V cxt (r)n(r)d 3 r 

t CW-/B«(M.A (la) 

where n(r) = (rJ(r), n^(r))(= (n(r),m(r))) labels the 
spin up and spin down charge density. 



E kin [n(r)} = ^<^a|T|^ a >, 



n(r)e xc [n{r)]d 3 r. 



(lb) 



(lc) 



Here T is an operator describing the kinetic energy and 
rest mass of the electrons 



rfr 

T= —a- V + mc 2 (/3-/ 4 ), 
i 



(2) 



where a and (3 are the usual relativistic matrices |4lj. 
U[n(r)) represents all two particle interactions including 
the Breit interaction. V CKt (r) is the external potential, 
B ext (r) is an external magnetic field. The density n(r) 
and the spin density m(r) are given by 



(3) 



II. THEORY 



m(r) = -^ B ^</'A( r )/ 3<T 4V'A(r), 



(4) 



A. The relativistic SIC-LSD formalism 

The SIC-LSD approximation |3ll l32j| is an ab-initio 
electronic structure scheme, that is capable of describ- 
ing localization phenomena in solids, [li* flfil Il7j In this 
scheme the spurious self-interaction of each occupied elec- 
tron state is subtracted from the conventional LSD ap- 



where 0*4 is the 4x4 matrix spin operator and A repre- 
sents the quantum numbers. In Eqs. and © below 
we have implied a representation in which spin is a good 
quantum number and the sums are over the occupied 
states. e;E C [rl(r)] is the exchange correlation energy of a 
gas of constant density and Eq. lfTc|) is the local spin 
density approximation. 
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If we minimise the functional (|la|) with respect to 
changes in the density and spin density we obtain a Dirac- 
like equation: 

ch 

(—a-V + mc^((3- h) + V cS (r) 
i 

+ MB/3cr 4 • B off (r))^ A (r) = e A ^A(r), (5a) 



where 



V ett (r) = V cxt (r) + 



n(r') 



3,,' 



47reo J |r — r'| 
6E^[n(r)} 



Sn(r) 



B cff (r) = B cxt (r) 



Sm(r) 



(5b) 
(5c) 



where n(r) = (n* (r),n^(r))(= (n(r), m(r))). The local 
spin density approximation discussed above provides a 
very successful description of a variety of properties of 
condensed matter, but suffers from a drawback because 
it contains self-interactions of the single particle charges. 
In an exact theory these spurious self-interactions would 
precisely cancel. In the LSD the cancellation is only ap- 
proximate and in materials where there arc well-localised 
electrons this can lead to significant errors. The SIC-LSD 
approach to this problem is to augment the LSD func- 
tional with an extra term that removes this deficiency. 



where 



E 



SIC-LSD 



E 



LSD 



E 



SIC 



(6a) 



£ SIC [{n 7 (r)}] = - $>K(r)] + £ x L c SD [n 7 (r)], (6b) 



where n 7 (r) = (rzl(r), ni ('))(= (%( r )> m 7( r ))) and 



U[nJr)\ = - 



1 e 2 



2 47T60 



n 7 (r)n 7 (r') , 3 , , 3 



' 7 ;, ' dVtfr, (6c) 
|r — r'| 



^ D [R 7 (r)] = j n 7 (r)e xc [n 7 (r)]d 3 r, (6d) 
where 7 runs over all orbitals that are Si-corrected, and 



ra 7 (r) = ^ 7 (r)^ T (r), 



m 7 (r) = -jjL-Bip\{Y)l3cTiip 1 {Y), 



(6e) 



(6f) 



by taking the functional derivative of E SIC LSD with re- 
spect to and we obtain 



, ch 



( — a-V + mc 2 {f3-h) + V ctt (r)+HBl3<T-B ctt {r) 

1 

+ y 7 SIC (r))^ 7 (r) = K?M r )> (7a) 

7' 



where the SIC potential is given by 



+ 



1 V 47r£ o 
«? x L c SD [n 7 (r)] 



n 7 (r') 



r — r'| 

^ D [n 7 (r)] 
5m 7 (r) 



(7b) 



The task of finding the single particle-like wavefunctions 
is now considerably more challenging than for the bare 
LSD because every state experiences a different potential. 
To maintain the orthogonality of the ip-y( r ) it is necessary 
to calculate the Lagrange multiplier matrix, A 77 ' . 

As written in Eqs. (10-13), £; SIC - LSD appears to be 
a functional of the set of occupied orbitals rather than 
of the total spin density only, like E hSr> . By a reformu- 
lation it may be shown Hi 111] that £ SIC - LSD can in 
fact be regarded as a functional of the total spin density 
only. The associated exchange-correlation energy func- 
tional E^. c [n(r)] is, however, only implicitly defined, 
for which reason the associated Kohn-Sham equations are 
rather impractical to exploit. For periodic solids the SIC- 
LSD approximation is a genuine extension of the LSD 
approximation in the sense that the self-interaction cor- 
rection is only finite for localized states, which means 
that if all valence states considered are Bloch-like single- 
particle states E SIC ~ Lsr> coincides with E LSD . Therefore, 
the LSD minimum is also a local minimum of ^sic-lsd 
In some cases another set of single-particle states may 
be found, not necessarily in Bloch form but, of course, 
equivalent to Bloch states, to provide a local minimum 
for £ SIC - LSD , For this to happen some states must exist 
which can benefit from the self-interaction term without 
losing too much band formation energy. This usually 
will be the case for rather well localized states like the 
3d states in transition metal oxides or the 4/ states in 
rare earth compounds. Thus, £; SIC - LSD i s a S pi n density 
functional, which may be used to describe localized as 
well as delocalized electron states. 

We have solved the SIC-LSD equations self- 
consistently for a periodic solid using the unified 
Hamiltonian approach described by Temmerman et at. 
[42j The equations have been solved on a periodic lattice 
using the relativistic LMTO method in the tight-binding 
representation. 



B. The relativistic spin-polarised LMTO method 



For the exchange-correlation term in the SIC energy we 
need to consider a fully spin-polarised electron. The cor- 
responding single particle-like wave equation is obtained 



In Section II C, u\' (r) will be a general notation for the 
unoccupied intermediate states in the second order time- 
dependent perturbation theory. In the case of a material 
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with translational periodicity U\> (r) will be a Bloch state 
u A ,(r) = ^ k (r), (8) 

for which 



^' k (r + R) = e lk - R ^' k (r), 



(9) 



where k is the wavevector defined to be in the first Bril- 
louin zone, j is the band index, and R is any Bravais 
lattice vector. In the LMTO method the Bloch wave 
functions may be expanded in several ways. |30| For the 
calculation of observables it is most convenient to make 
an expansion in terms of the single-site solutions of the 
radial Dirac equation and their energy derivatives. For 
the relativistic spin-polarised case this has been achieved 
by Ebert [43l HI ] and it is this method that we employ. 
The Bloch state in this representation is written as 



V> jk (r) 



JVtypo N t 

t=l i=l A 



<AWr-?f )]• 



(10a) 



Here 



where the g^}^ 1 {E 
of the radial Dirac equation for a spin-polarised system 
and rj^ = t — t^. Details of the solution are given in by 

Strange et ai, |45| and </> i ,tA( r i*'') i s its energy derivative. 
These satisfy 



^\if^^rf)^(^))' (1 ° b) 



r.-*^) and f^™ 1 '(E,r^) are solutions 



HA 



jtA) 



1 



W|^ tA )=0, (11) 



where the subscript v corresponds to the energy e v about 
which the muffin-tin orbitals of Eq. IjlObfl are expanded, 
and the normalization integrals have been done within 
the atomic sphere 5 1 *. The single particle functions 

</wa(i") and </>„ tA (r) are evaluated at energy e v . In this 
relativistic formulation A = (nrrij) labels the bound- 
ary condition for the independent single-site solution 



<?W(r- 



.(«)x 



at r. 



(t) 



N 



type 



of the Dirac equation about the basis atom 
is the number of different types of atom 



in the unit cell. Nt is the number of equivalent atoms 
of type t. The coefficients A^ iA and -B^A are written 
in terms of the LMTO structure constants and poten- 
tial parameters, and are completely determined by a self- 
consistent LMTO calculation of the electronic structure. 
[3fj| Key observables are then given in terms of these 
quantities. In particular the spin moment is 



ms 



(Tk, 



eJ k <e F 



(12a) 



where 



t,i A A' 



+ B tiA*- B MA'(^tAl/ 3o "4 Z |^iA'! 
+ B tiA A tiA>(<t>utA\Pa±z\<l>vtA> 



(12b) 



with ep being the Fermi energy and e jk is the Bloch state 
eigenenergy. The orbital moment is 



m L 



E 



(13a) 



where 



t.i A A' 
?ik* R jk 



+ Ka*Ka> (<f>utA\Plz\<t>utA') 
+ Bi^A%(^ tA \(3l z \M 



(13b) 



In all our calculations the B-field is along the z-axis which 
therefore acts as an axis of quantization. 



C. The x-ray scattering cross section 

In this Section we will outline the formal first- 
principles theory of magnetic x-ray scattering for materi- 
als with translational periodicity. The theory is based on 
the fully relativistic spin-polarised SIC-LMTO method 
in conjunction with 2nd order time-dependent pertur- 
bation theory. To simplify the presentation a straight- 
forward canonical perturbation theory [41| is presented 
rather than a more sophisticated diagrammatic method. 

M 



1. Basic theory of x-ray scattering 

The theory of x-ray scattering is based on the second 
order golden rule for the transition probability per unit 
time: 



"if 



2tt 



(WL 



E 



(f\H! nt \I)(I\H!Ji) 



E, 



-Ei 
xS(E f 



Ei) (W) 



where \i), \I) and |/) are the initial, intermediate and fi- 
nal states of the electron-photon system. Ei, Ej, and 
Ef are the corresponding energies. H[ nt is the time- 
independent part of the photon-electron interaction oper- 
ator. The formalism to reduce this general expression to 
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single-electron-like form has been published previously. 
|lOl | Therefore we will not repeat the details here, but 
only the equations that are key to the present implemen- 
tation. 

In relativistic quantum theory it is the second term 
in Eq. I|15|) that is entirely responsible for scattering as 
it is second order in the vector potential. It is conve- 
nient to divide this term into four components. To see 

I 



where u\(r) and i>a(i") are positive-energy electron and 
positron eigenstates of the Dirac Hamiltonian for the 
crystal and form a complete orthonormal set of four- 
component basis functions in the Dirac space. The quan- 
tum state label A can then be related by symmetry ar- 
guments to A. In Eq. (|15|l term [1] represents scattering 
with no photons and positive energy electrons only in the 
intermediate state, term [2] is when there are two pho- 
tons and positive energy electrons only in the intermedi- 
ate state, term [3] is for no photons and when negative- 
energy electrons exist in the intermediate state and term 
[4] is for when two photons and negative-energy electrons 
exist in the intermediate state. We may recall that within 
the golden rule based Thomson scattering formalism the 
negative-energy related state terms have the wrong sign. 
Therefore amplitudes [3] and [4] have been non-rigorously 
corrected by multiplying them by —1. The positive en- 
ergy one-electron states are subject to the constraint that 
(a < £f and e\i > ep. The relativistic photon-electron 
interaction vertex is 

/ fir 2 \ 1/2 

A qA (r) = - e — — a.eW(q)e*- (16) 
\2Ve LuJ 

where e = — |e|, and q, A (q', A') represent the wavevector 
and polarisation of the incident (outgoing) photon, and 
e( A ' (q) is the polarization vector for the x ray propagat- 
ing in the direction of q. The a = (a x , a y , a z ) are the 



this note that there are just two types of intermediate 
state \I), those containing no photons and those contain- 
ing two photons. We can also divide up the scattering 
amplitude according to whether or not the intermediate 
states contain excitations from the 'negative-energy sea 
of electrons', i.e. the creation of electron-positron pairs. 
It can be shown that the x-ray scattering amplitude in 
the case of elastic scattering can be written as [ljj, |46J 



[1] 
[2] 
[3] 

W, (15) 



usual relativistic matrices in the standard representation. 
In Eq. (|15fl the last two terms are neglected. The justi- 
fication for this is twofold. Firstly, in the energy range 
of interest hw « 2mc 2 these two terms have no reso- 
nance, and so will only make a contribution to the cross 
section that is slowly varying. This is to be compared 
with the resonant behaviour of the first term. Secondly, 
in Thomson scattering, where the negative energy states 
play a key role, all the electron states are extended. In 
a crystalline environment the negative energy states are 
largely extended while the states close to the Fermi en- 
ergy are more localised, so one would expect the matrix 
elements to be smaller. For further details see Section II 
C of Ref. ^3 Henceforth the first term in Eq. l(T5|) will 
be referred to as the resonant term and the second as the 
non-resonant term. 

In elastic scattering of x rays u A (r) is an atomic-like 
core state localised at a lattice site. Although it is lo- 
calised it is still an electron state of the crystal Hamilto- 
nian. It is given by 

UAjri) "tWS7(^)x^(f n )J' (17) 

where g^™ 1 (r n ) and /^"i™ 3 ( r «) are solutions of the ra- 
dial spin-polarised Dirac equation |45| at the site n and 



/qA;q A M = ^H + ^H + Q A ,H + Q,H 



y, (J\H! nt \I)(I\H! mt \i) 

/,CA>0 



y 

E * E ' ,,tr<o Ei - Ei 

J d 3 ru{(r)xl, x , (r)u A , (r) J d 3 r'u{, (r')A qA (r> A (rQ 
£A — eA' + fuv 

/d 3 TO i(r)A qA (r)u A ,(r)/d 3 r^ A ,(rOA q<v (rQ MA (r') 

£A -CA> -tUV 

/rf 3 r4(r)At <A< (r)^(r)/d 3 r-u A ,(rOA qA (rO^(rO 

£A — £A' + 

J d 3 r4(r)A qA (r)u A , (r) J d*r'u\, (r')A q , A , (r'^r') 

£A - £A' - t^> 

I 



(f\H! nt \I)(I\H! nt \i) 



E 

AA' 

E 

AA' 

E 

AA' 

E 

AA' 
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Xk 3 (r) are the usual spin-angular functions with angu- 
lar momentum related quantum numbers A = (nrrij). 
[ill ItH EsT ] As in Eq. I|10b(l the sum over k' u runs over 
k' = n n and n' n = —Kn - 1 only. 



U. Evaluation of the cross section 

The physical observable measured in MXRS experi- 
ments is the elastic differential cross section for scatter- 
ing. This is given by (see Sec. II E of Ref. ITfl) 



2, ,2 



da 



dtt (2n) 2 h 2 c z 



|/qA;q'A'H| , 



(18) 



where the symbols have their usual meanings, and we 
need to calculate the first two terms of Eq. Q15[l. i.e. 

j+(pos) anc | j-(pos)^ 

When implementing Eq. i|15fl for a perfect, translation- 
ally periodic multi-atom per unit cell crystal we use the 
following coordinate transformations 



r = R/, + r/, = R,f 3 + R/ + v h , 
r' = Rj. + r j. = Rf + R, 7 + r. 7j 



(19a) 
(19b) 



where r|°' > and denote the i th and j th basis atoms, 



respectively, in the th unit cell, and R 7 and Rj are 
Bravais lattice vectors. 

Furthermore, we use the substitutions 



E 

A 

E 



EE E . 

W " A(A„) 

E< 



and 



'3,..' 



•/ 3 



(19c) 
(19d) 



(19e) 
(19f) 



where N, I and J stand for the label of unit cells, n, i, 
j stand for the label of basis atoms, and A(N n ) is the 
initial core state label for an atom at site Re- 
using Eq. I|19fl and Eq. © in connection with term 
[1] of Eq. l(TB|) the resonant part of the positive-energy 
scattering amplitude for a perfect crystal can be written 
as 



f +( P os) , N 

/qA;q'A' V W ) 



EE EE 

N n A(N n ) jk 



e A(A„) 



EE / ^ u A(A„)(R/ 1 +r/J^ q , A ,(R /! +r /i )V jk (R/ 1 +r/,) 
EE /, rfy 7j ^ k ^R/ J +r' Jj )^qA(R./ J +r: 7j ) UA(Wn) (R j3 +r' Jj ), 

7 „- JS 1 



(20) 



where the sums are restricted such that £a(./v„) < £f and 
e jk > e F . 

We approximate Eq. (|20(l in a similar way as we did 
earlier in our R-SP-GF-MS method based MXRS theory 
(see Section II B of Arola et al. 1 Oj ) . Because the core 
states ma(tv„) participating to the x-ray scattering (XS) 
are well-localized around site Rat„ , the dominant contri- 
bution to XS in Eq. H2()|) becomes from the Ii — N n (i.e. 
I = N, i = n) and Jj — N n (i.e. J = N, j — n) terms. 
From the physical viewpoint, this refers to the situation 
where in the anomalous scattering process of x rays a 
core electron will be annihilated and created at the same 
atomic site (site-diagonal scattering). 

Furthermore, we note that in the perfect crystal case 
the following properties can be used: 1) Rw„ = Rn°' ) + 
Rat; 2) electronic coordinate rjv„ can be replaced by r„; 



3) V ik (RT + Ra + r„) = e* k - R «^' k (Rr + r„), i.e. 
Bloch's theorem for intermediate states; and 4) the core 
state label A(N n ) = A n , i.e. is unit cell independent. If 
we also use the explicit form of the photon-electron inter- 
action vertex of Eq. I|16fl then we end up to the following 
expression for the /+(p os ) scattering amplitude: 



«h = E EEE^ (q '~ q)R 

AT I jk n A„ 



(0) 



4';)^ k (q'A')mi"j +jk *(qA) 



» 



a -i(q'-q)Rjv 



hhJ 



(21a) 



7 



where the resonant matrix elements are defined as 



An 



(qA) = J s J 3 r n u^Hr n )xl x (r n )^ k (R^+r n ), 

(21b) 

where S n refers to the n th atomic sphere within the unit 
cell. 

In Eq. (|21|l we notice that 



£- e -i(q'-q).RAT = ^ £ ^_ q , K , 



JV 



K 



and 



(22a) 



(22b) 



where N ce \\ s is the number of unit cells in the crystal, K is 
a reciprocal lattice vector, and <5 q '_ qi K is the Kronecker 
6 function. 

As the last step, we decompose the general basis atom 
label n in Eq. (|21|l into type t it = 1, . . . , N ty p B ) and basis 
atom label i (i = 1, ... ,N t ), i.e. n = (ti). Consequently, 
this introduces the following notational changes in Eq. 

n t=l i=l 

(*) _> „(*) v. . „(*) 



1/ 



u w ,r fe -> r 



(23) 



Implementing these notations along with Eq. i|22|) in 
Eq. (|21|1 . leads to the final expression for the resonant 
part of the scattering amplitude in Bragg diffraction 
which is 

= f£$$A»)Nc*. y E6 Q K, (24a) 

K 

where the th unit cell contribution to the scattering am- 
plitude is 

.„ (t) mf )+jk (q'A')™l tl)+jk *(qA) .. 



e A< £j ' k + ^ + ir A*/ 2 



(24b) 



where Q = q' — q, and the matrix elements m A "^(qA) 
are given by Eq. i|21bl) with the new notations of Eq. (|23|1 . 

The added phenomenological parameter represents 
the natural width of the intermediate states created by 
the core hole state |A t ) at the t type basis atom. 

Similarly, starting from term [2] of Eq. (|15|l , it can be 
shown that the expression for the nonresonant part of the 



scattering amplitude in Bragg diffraction can be written 
as 

f&lM = /oTSU^cellsE^QK, (25a) 

K 

where the th unit cell contribution to the Bragg scatter- 
ing amplitude is 



/0:qA°q'A'M - E 



V 



,.>_,:; / d * k E EE 



xe 



W mi t f- jk (qA)mi*; ) -^(q'A') 



e £ - eJk - fiw 



(25b) 



where the non-resonant matrix element is defined as 



(25c) 



The total amplitude in Bragg diffraction can then be cal- 
culated from 



K 



(26) 

where fo^ poB ^ and / ' pos ' ) amplitudes are given in Eqs. 



24b(l and l|25b|) . respectively. 



D. Matrix elements 

In this section we present the derivation of computable 
expressions for the matrix elements ro A +jk (qA) and 
m («)-jl*( qA ) in the f ram ework of the R-SP-SIC-LMTO 
electronic structure method. 

Using the expansion of Eq. Hl()a(l for the SIC-LDA 
Bloch state ip 3 (r), and noticing that the independent 
single-site solution </vtA(r) of the Dirac equation vanishes 
outside the atomic sphere at site (ti), i.e. 0„tA( r ) = 
for r > S* (cf. Ref. l49l pp. 120-1), then the resonant 
matrix element m A tI '' +: ' k (qA) can be written as 



in 



(H)+jk 
A t 



(qA) = EK'a^I^aI^a) 



tA, 



(27a) 



where (f\X^ x \g) is defined as 



(f\xl x \g)= f d 3 r/t(r)X qA (r) 9 (r), (27b) 



where / = u£ and # = <j> vtA or <?W. 
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Similarly, it can be shown that the nonresonant matrix 
elements can be written as 



(28) 



where (f\Xq\\g) quantities can be calculated by doing 
the replacement X^ x — > Xq\ in Eq. J27EJ). 

Finally, we mention few practical points about the im- 
plementation of the matrix elements of Eqs. (|27[1 and 
(|28[1 . We will derive below numerically tractable approx- 
imations for these matrix elements due to the electric 
dipole (El) or magnetic dipole and electric quadrupole 
(Ml + E2) contributions to the photon-electron interac- 
tion vertex Xq\(r). 



in terms of the radial and angular integrals; Rws is the 
Wiegner-Seitz radius and the angular integrals are de- 
fined by [13 



A 



Km, 7 \K'm-i 



X™ 3t (r)<x 



(q) X 7(f)da (30b) 



A numerically tractable expression for 
&\Xa) )\<j) u th) in Eq. (|29af) can then be written 



(t)rrij 



immediately by doing the replacements f^, K 
9k'k 1 on tne right side of Eq 



and &* 



1. Matrix elements in electric dipole approximation 



In the electric dipole approximation (El) [e lq r w 1 in 



Eq. 1)16(1]. the resonant matrix element 



l At 



(qA) Of 



Eq. H27[l can be written as 

maf?+* (qA ) = E[ A ua(4!\H>\M 

A 

+ B^(uf t \Xa[ x \^ tA )\ , (29a) 

where 



Xa,q X (r) 



-ec 



(— 

\2Ve io q 



1/2 



a-eW(q). (29b) 



Using the core state expansion, Eq. i|17|). and the 
expansions l|10b|) for <p ut A( r ) an d the analogous ex- 
pansion 4> v t a(y)i respectively, it can be shown that 
(u At \XaL x \(f) v t\) in Eq. I|29a|l can be written as 



( u A t l Xa qAl<?W) = ~iec 



qAlT ' ' ~\2Ve uJ q 



1/2 



E 

k 4 'k' 



1 ws 





(-A) 



x A K ~" J (a) 

Kt'Vflj t \ — K'mj 



x A 



o 

l (_A ) - (q)l 



drr fJ< Kt W K ' K 3 (tu,r) 



(30a) 



Similarly, making a replacement Xq\ —> Xa,q\ in Eq. 
and using the property Ia q j(r) = Xa^_ x (v) for 
circularly polarized light, we can immediately show that 
the nonresonant matrix elements can be computed from 
the resonant ones in the El approximation (for further 
details, see Sec. II D of Rcf. lit; as 



(ti)-j'k 



(qA) = 



fti)+jk 



(q-A). 



(31) 



If the photon propagates along the direction of magneti- 
sation (the z-axis) then the unit polarisation vectors for 
left (LCP) and right (RCP) circularly polarised light are 
e(+)(z) = (l,i,0)/ v / 2ande(-)(z) = (1, -i, 0)/\/2 respec- 
tively. To obtain the polarisation vectors for propagation 
directions away from the z-axis rotation matrices are ap- 
plied to these vectors. Using the well-known orthonor- 
mality properties of the spherical harmonics the angular 
integrals of Eq. (|30bjl can be written 
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+ f2i{6 q ,<j> q ,\)C(l^j;m j + ^,-i)C(7'i/;m^ - ^, ^M'^m^-i 

+ /22(^,^,A)C(/ii;m j + I,-i)C(l'i/;m5 + 5, -^^j. (32) 



The angular factors fij(9, 4>, A) are determined by the di- 
rection of propagation and the photon polarisation. They 
are discussed in detail by Arola et al. 01 I* 1 the case 



where the direction of q is described by a rotation around 
the j/-axis of q followed by a rotation about the z-axis 
of <j) q , in the active interpretation they are given by 



/n( 

/2l( 

fn( 

/2l( 



'q,Yq, 
\Aqi- 

\Aqi- 
^q,4>q,- 



-) fl2(0q,<Pq,+) 

-) fa(8 9 ,<t> 9 ,+) 

-) fl2(6q,(j> q ,-) 

-) fa(fi 9 ,<t> 9 ,-) 



If — sin 9 (cos # g 

71 v(cos6» g - l)exp(i^) 

1 f — sin g (cos ^ 



+ 1) exp(-i0q)\ 
sin # g y 



sin 

- 1) exp(-i0 9 
sin 9 q 



(33) 
(34) 



for positive and negative helicity x rays, respectively. The 
angular matrix elements of Eq. I)30bjl together with the 
symmetry of the A J tiA and B J UA coefficients determines 
the selection rules in the electric dipole approximation. 

It is important to note that the selection rules, derived 
originally for x-ray scattering in the framework of the 
Green's function multiple scattering electronic structure 
theory, 0] can be applied as such only to each term of 
Eq. (I29a.fl separately with angular momentum -like quan- 
tum numbers of the core state (A t ) and single-site valence 
orbital (A). [5(| The El selection rules then become 
I — It = ±1 for RCP and LCP radiation in any propa- 



gation direction, while mj 



m 



0, ±1, depending on 



the polarization state as well as on the propagation di- 
rection of the photon. It is also noticeable that the selec- 
tion rules in the case of matrix elements ma A ^ U (qA) 
are slightly different from the case of matrix elements 



ma 



At 



(qA) with respect to the azimuthal rrij quan- 
tum number, because Eq. (|30afl contains angular matrix 
elements of the form A(~ x \ while the corresponding ex- 
pression for ma^~"' k (qA) contains with an opposite 
polarization state index. [T(i| 

Derivation of the selection rules for the matrix elements 
ma^ t ^ +: ' k (qA) or ma A ^ ~"' k (qA) would be possible only 
for k points of high symmetry whose irreducible double 
point group representations and the angular momentum 
A (nrrij) decomposition for their symmetrized wavefunc- 
tions are known. However, we apply numerical rather 
than group theoretical procedure to determine the selec- 
tion rule properties of the abovementioned matrix ele- 



ments. 



2. Matrix elements due to magnetic dipole 
and electric quadrupole correction 

We derive below an expression for the combined mag- 
netic dipole and electric quadrupole (M1+E2) correction 
to the electric dipole approximation (El) of the matrix 
elements of Eqs. (|27aH and (|28fl . If we now approximate 



1-Mq-r in Eq. (|16fl for Xq\(r), then the term iqr 
is responsible for the (Ml + E2) corrections to the elec- 
tric dipole approximated matrix elements ma A ^ (qA) 



and ma 



A t 



and mb 



(ti)-jk 

At 



(qA), which we denote as mb 
(qA), respectively. 



"A t 
At 



(qA) 



It is then a straightforward matter to show that the 



matrix element 



(qA), related to the resonant 



part of the scattering amplitude, can be written as 



mb 



A t 



(qA) 



E 

A 



[A^(uf t \Xb[ x \^ tA ) 



B&{uV>\Xbl x \M , (35a) 



where 



A& qA (r) 



-ec 



2Ve uj q 



1/2 



(q)iq • 



(35b) 



Using again the angular momentum expansions of the 
,tA, and <j) u tK functions, as we did in the 



core state u 



(*) 

At' 
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derivation of the ma^ + ^ (qX) expression, we get for the 
first term of Eq. as 



{ufAXb^M = -ec 



2Ve uj q 



1/2 



i?, 





(-A) 



x H y ~^' (a) 



x B 



o 



t' rrij t \K,'m,j \1/ J ' 



(36a) 



where the angular integrals are defined by [Ti 



d(A) /^n 



m,-f/ 



(f)(T.e( A )(q)q.f X 7(f)^, 



(36b) 
where |q| = |f | = 1. 

A similar expression can be worked out for the second 

term of Eq. Ij35a(l by doing the replacements f^)™ 1 — > 

-> 9^./ ™ 3 on the right side of Eq. (|36all . 
the nonresonant matrix elements 



Jk'k ancl 9 K ' K 
By Eq. QZa>, 



to6 



(ti)-jk 

At 



(qA), due to the (M1 + .E2) correction, can then 



be worked out by making the replacement Xb^ x — ► A6 q A 
inEq. I|35a(l . By noticing that X6 q },(r) = — A& q _ A (r), we 

£72) matrix elements 



can express the nonresonant (Ml 
in terms of the resonant ones as 



mb 



(ti)-jk 

At 



(qA) 



, (ti)+jk, 

-mb\ t > (q 



A). 



(37) 



The angular matrix elements B' (a) of Eq. 136bll 

j > j 

can be written as a sum of twelve terms (see Eq. (26) of 
Ref . IT(ih . Consequently, the selection rules of the (Ml + 
E2) contribution to the x-ray scattering are essentially 
more complicated than in the El case. 

As guided by the El case above, we can derive the 
(A/1 + E2) selection rules for each term of Eq. Ij35a(l 
with angular momentum -like quantum numbers of the 
core state (A t ) and single-site valence orbital (A). The 
resulting selection rules are then I — It = 0,±1,±2 
with the restriction that s — > p and p — > s be forbid- 
den transitions, and for the azimuthal quantum number 
TOj — rrij. = 0,±1,±2, depending on the direction and 
polarization of the photon. ^(J 



III. RESULTS 

In this section we discuss a series of calculations to il- 
lustrate the relativistic MXRS theory we have developed 
within the SIC-LSD method for ordered magnetic crys- 
tals, and to demonstrate explicitly what information is 
contained in the x-ray scattering cross section. For this 



en 

CD 
_Q 

>N 



O 

o 



-4.698 




-4.702 



Wigner Seitz Radius (a.u.) 

FIG. 1: The calculated SIC-LSD total energy of fee 
praseodymium as a function of Wigner-Seitz radius. The elec- 
tronic configuration corresponding to the Hund's rule ground 
state was used for these calculations. The theoretical predic- 
tion of the Wigner-Seitz radius is 3.82 a.u. Experimentally 
(see text) the value is 3.818 a.u, |3(J or 3.827 a.u. (5l| or 
3.793 a.u. El. 



we have chosen to examine fee praseodymium for a de- 
tailed analysis of the theory. The reasons for this choice 
are as follows: (i) Praseodymium contains two localized 
/-electrons. Therefore, it is the simplest /-electron mate- 
rial for which we can to a large extent alter both the spin 
and orbital contributions to the magnetic moment by se- 
lectively choosing beforehand for which electrons we ap- 
ply the SIC correction, (ii) Being ferromagnetic and fee it 
has only one atom per primitive cell and is therefore com- 
putationally efficient to work with. While the fee struc- 
ture is not the observed ground state of Pr, it has been 
fabricated with this structure at high temperatures and 
pressures, (iii) Using nonrelativistic SIC-LSD we have 
obtained good agreement with experiment for the valence 
and equilibrium lattice constant of praseodymium, (iv) 
Preliminary calculations indicate that for the rare earth 
Miv and My edges the MXRS spectra are, to first order, 
independent of crystal structure, so the results we obtain 
may be provisionally compared with experiment. 



A. Ground state properties 

We have performed a self-consistent fully relativis- 
tic SIC-LSD calculation of the electronic structure of 
praseodymium at a series of lattice constants on the fee 
structure and found a minimum in the total energy as 
shown in Fig. ^ where the results are presented in terms 
of the Wigner-Seitz radius. There is are variety of differ- 
ent methods for obtaining the experimental lattice con- 
stant. Firstly we can use the Wigner-Seitz radius that 
corresponds to the same volume per atom on the fee lat- 
tice as is found in the naturally occuring dhep crystal 
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TABLE I: This table displays the /-states selected for self- 
interaction correction and the self-consistently calculated spin 
and orbital magnetic moments of those states. The first col- 
umn simply labels different configurations of localized states, 
the second column gives the mi and m s quantum numbers 
of the states from which self-interactions have been removed. 
Columns 3 and 4 are the calculated spin and orbital contribu- 
tion to the total magnetic moment from the self-interaction 
corrected /-electrons shown in column 2. Columns 5 and 6 are 
the calculated total spin and orbital contribution to the mag- 
netic moment from all electrons in fee Pr. Note that the spin 
moment is fairly constant for all the selected configurations. 



0.6 



N 


(mi,m s )i 


(m(,m s )2 


1VL (SIC) 


M ; (SIC) 


M L (t) 


M s (t) 


1 


(-3,1/2), 


(-2,1/2) 


-4.96 


+1.98 


-4.79 


2.42 


2 


(-3,1/2), 


(-1,1/2) 


-3.95 


+2.00 


-3.91 


2.43 


3 


(-3,1/2), 


(0,1/2) 


-2.99 


+1.99 


-2.97 


2.43 


4 


(-3,1/2), 


(1,1/2) 


-1.96 


+1.96 


-1.97 


2.48 


5 


(-3,1/2), 


(2,1/2) 


-0.98 


+1.99 


-1.03 


2.48 


6 


(-3,1/2), 


(3,1/2) 


-0.005 


+2.00 


-0.05 


2.52 


7 


(3,1/2), 


(-2,1/2) 


1.00 


+1.99 


0.89 


2.48 


8 


(3,1/2), 


(-1,1/2) 


2.02 


+1.97 


1.84 


2.50 


9 


(3,1/2), 


(0,1/2) 


3.00 


+1.98 


2.75 


2.49 


10 


(3,1/2), 


(1,1/2) 


4.00 


+1.97 


3.80 


2.52 


11 


(3,1/2), 


(2,1/2) 


4.99 


+1.99 


4.69 


2.54 
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Orbital Moment (Bohr Magnetons) 

FIG. 2: The calculated total energy per atom of fee 
praseodymium relative to the ground state energy as a func- 
tion of orbital moment in Bohr magnetons for the states 
shown in table I. All magnetic moments were self-consistently 
determined and the spin moment was approximately constant 
for all the configurations shown. If an antiparallel arrange- 
ment of spins was selected the energies were considerably 
higher. 



structure [3(j • This gives i?ws = 3.818 a.u. Secondly 
we can take the room temperature value which is ob- 
tained experimentally from flakes of Pr by quenching in 
an arc furnace. This yields i?ws = 3.827 a.u. and we can 
take the value reported by Kutznetsov [E^, i?ws = 3.793 
a.u. which was measured on samples at 575 K. Clearly 
our calculated value of i?ws = 3.82 a.u. is in excellent 
agreement with these values. Following earlier work by 
Myron and Liu |53| Soderlind performed a comprehen- 
sive first-principles study of the electronic structure of 
Pr using the full potential LMTO method which shows 
that the f ee p hase is stable at pressures between 60 and 
165 kbar yy|. Calculations employing the SIC within a 
non-relativistic framework have been performed by Tem- 
merman et al. j3|| and by Svane et al. [5^ 

Within the SIC-LSD method we can choose which elec- 
tron states to correct for self-interaction. As the effect of 
the SIC is to localize the states this effectively determines 
which two of the 14 possible / states are occupied in 
trivalent praseodymium. All non-SI-corrected electrons 
are described using the standard local spin density ap- 
proximation via the unified Hamiltonian describing both 
localized and itinerant electrons. By trying all possible 
configurations and determining which arrangement of / 
electrons has the lowest total energy we can determine 
the ground state of praseodymium. It should be pointed 
out that this interpretation is rather distinct from the 
standard model of the rare earth magnetism where the 
Hund's rule ground state can be thought of as a linear 
combination of possible 4/ 2 states. In our model the 
exchange field is automatically included and this yields 
a Zeeman-like splitting of the 4/ 2 states and gives us a 



unique ground state. In Table I we show a selection of 
possible states occupied by the two electrons with their 
self-consistently evaluated spin and orbital magnetic mo- 
ments. In Fig. 0we display the calculated total energy of 
these states against orbital moment. Note that we have 
chosen the spin moments parallel for all the states shown. 
For the antiparallel arrangement of electron spins the en- 
ergy is significantly higher. It is clear that there is an ap- 
proximately linear relationship between the total energy 
and orbital moment. For all the points on this figure 
the orbital as well as spin moments are computed self- 
consistently including the relaxation of the core states. 
The spin moments for all configurations of / electrons 
are found to be approximately the same, always being 
within 0.06 /j,b of 2.48 /ib (See Table I) in fair agreement 
with the result of Soderlind 541]. There is also a small 
increase in the magnitude of the (positive) spin moment 
as the orbital moment increases from its most negative 
to its most positive values. This is due to the increas- 
ing effective field felt by the valence electrons. There 
is a slight variation in the spin moment values because 
the small hybridization of the non-SIC corrected / elec- 
trons with the 5d — 6s conduction band is dependent 
on the orbital character of the occupied states. Fig. [21 
is consistent with the Hund's rules. The lowest energy 
state has /-spins parallel to each other in agreement with 
Hund's first rule. The total spin moment is 2.42 /ib of 
which the two localized / electrons contribute 1.98 fiB, 
and the remainder comes from spin-polarization in the 
valence bands. The z component of the orbital mag- 
netic moment is -4.79 which is composed of -4.97 fj,B 
from the localized /-electrons and 0.18 /xb from the va- 
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lence electrons. Note that the valence contribution to the 
orbital moment is parallel to the spin moment and an- 
tiparallel to the localised orbital moment. These numbers 
are fully consistent with Hund's second rule. Further- 
more the /-shell is less than half full and the spin and 
orbital moments are found to be antiparallel in the low- 
est energy state, consistent with Hund's third rule. The 
fact that we can reproduce the expected lattice constant 
and /-electron configuration suggests very strongly that 
the electronic structure calculated using the relativistic 
SIC-LSD method describes the ground state properties 
of fee praseodymium well. A detailed discussion of the 
electronic structure of the rare earth metals, calculated 
using the relativistic SIC-LSD method, will be published 
elsewhere [Hfij 

B. X-Ray scattering cross sections 

We have performed calculations of the x-ray scatter- 
ing cross section at the Mjy and My absorption edges 
of Pr for all the /-electronic configurations shown in Ta- 
ble I. These were evaluated with an arbitrary value of 
r=l eV which is smaller than one would expect experi- 
mentally, but as the purpose of this section is to inves- 
tigate the capability of the theory only rather than to 
make a strict comparison with experiment, it does not 
pose a problem. The effect of increasing T is simply to 
broaden and smooth out the calculated curve. A selec- 
tion of the results is shown in Figs. 01 and 0] In each 
figure the cross section for left-handedly circularly po- 
larized (LCP) photons, i.e. with positive helicity, and 
right-handedly circularly polarized (RCP) photons, i.e. 
with negative helicity, arc shown. The geometrical setup 
of these calculations assumes that the propagation direc- 
tions of the incident and outgoing photons are parallel 
to the exchange field, i.e. in our case parallel to the spin 
magnetic moment. Figures and 0] show the cross sec- 
tion at the Miv and My edges respectively as the SIC 
configuration is changed systematically such that the z 
component of the orbital moment varies from negative 
to positive values while at the same time the calculation 
shows that the spin moments remain nearly constant in 
magnitude and parallel to the exchange field. As the or- 
bital moment increases we see that the cross section at 
the Miv edge changes only slightly for LCP x rays while 
for RCP x rays it changes dramatically. At the most neg- 
ative orbital moment the Miv RCP cross section is very 
small, being completely overshadowed by the LCP peak. 
At the other end of the scale where the orbital moment 
is most positive the cross section for Miv RCP x rays is 



considerably larger than that for LCP x rays. It should 
also be noted that the cross section peak for RCP x rays 
is 1 — 2 eV lower in energy than the peak for LCP x rays. 

When the resonant scattering (hut w ep — £A t ) is close 
to the My edge, it is the RCP cross section that remains 
approximately constant with changing orbital moment, 
although a significant shoulder does appear on the low 
energy side of the curve as the orbital moment increases. 
The LCP peak decreases dramatically with increasing or- 
bital moment. At the My edge, peaks from RCP and 
LCP x-ray scattering are again separated by 1 — 2 eV, 
but the ordering of the peaks is reversed from the case of 
the Miv edge scattering. 

Figures and 0] indicate that the Miv and My cross 
sections are directly related to the orbital moment of the 
constituent atoms, although they do not indicate the di- 
rect proportionality between magnetic moment and scat- 
tering cross section suggested by Blume For example, 
the My edge cross section for LCP photons hardly varies 
in the upper two pictures in Fig. 0] despite a change of 
nearly 2fiB in the orbital moment. To clarify this point 
further, we show in Fig. [S]the cross section at the Miv 
and My edges for SIC configurations that produce an 
orbital moment close to zero with the spins of the two 
occupied / states parallel. While neither the spin nor 
the orbital moment change significantly, the cross sec- 
tion certainly does. At the My edge the negative helicity 
curve is approximately constant while the positive helic- 
ity curve alters dramatically. On the other hand, at the 
Miv edge it is the positive helicity curve that is approxi- 
mately constant while the negative one shows significant 
variation. This figure implies that the resonant x-ray 
scattering does not measure the total orbital moment, 
but is a measure of the orbital angular momentum of the 
individual one-electron states. 

The important message of Figs. OHSlis that the scatter- 
ing cross section is not directly proportional to the total 
orbital moment of the material. However, both the spin 
and orbital moment have a strong influence on the size 
of the cross section peaks. 



IV. DISCUSSION 

The standard theory of x-ray magnetic scattering is 
based on the work of Blume Q] . He derived an equation 
for the nonresonant x-ray scattering cross section using 
a nonrclativistic approach with relativistic correction to 
order 1/c 2 . The resulting expression for the cross section, 
using his notation, is 
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FIG. 3: The scattering cross section at the Mrv edge for praseodymium for electron configurations 1,3,5,7,9, and 11 from 
Table I. Each figure is for a different pair of localized /-electrons. The calculated total orbital moment in Bohr magnetons is 
shown in the top left of each figure. The full curve is the cross section for x rays with positive helicity and the dashed curve 
is that for negative helicity x rays. A general trend of increasing magnitude of the cross section for negative helicity incident 
photons as the orbital moment increases from negative to positive is clearly observable in these curves. The positive helicity 
curve remains approximately constant with increasing orbital moment. 
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where E a and E b are the energies of the initial and fi- 
nal many-electron states, |o) and respectively. K = 
k k', where k and k' are the wave vectors of the in- 
coming and scattered photons, respectively, and rj, Pj, 
and Sj are electron (rigorously density functional state) 
coordinate, momentum, and spin operators. 

A = e x e (39) 

and 

B = e' x e - (k'xe')x(kx e) 

- (k x e)(e' ■ k) + (£:' x e')(e ■ k') (40) 

depend only on the direction and polarization of the in- 
cident and emitted photons. The first term in Eq. Ij38(l is 
the Thomson term, responsible for the charge scattering. 



The term containing A depends on the orbital momen- 
tum and the term containing B depends on the elec- 
tron spin. This expression clearly shows that there are 
three distinct contributions to the magnetic scattering 
cross section, one from the orbital moment, the second 
from the spin moment, and the third from the interfer- 
ence term between the spin- and orbital moment. We 
also note the obvious point that if the orbital and spin 
moments of the individual electrons sum to zero, then 
the magnetic scattering vanishes. Most interestingly, Eq. 
(I38JI implies that, with a suitable choice of the photon 
energy, geometry and photon polarization it is possible 
to separate contributions to the cross section from the or- 
bital and spin moments. However, this expression is not 
directly applicable in our resonant magnetic scattering 
studies because its derivation involves an approximation 
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FIG. 4: The scattering cross section at the My edge for praseodymium for electron configurations 1, 3, 5, 7, 9 and 11 from Table 
I. Each figure is for a different pair of localized /-electrons. The calculated total orbital moment in Bohr magnetons is shown 
in the top left of each figure. The full curve is the cross section for x rays with positive helicity and the dashed curve is that for 
negative helicity x rays. A general trend of decreasing magnitude of the cross section for positive helicity incident photons with 
increasing orbital moment is clearly observable in these figures. The negative helicity curve remains fairly constant in magnitude 
with increasing orbital moment although the feature on the low energy side of the peak does become more pronounced. 



which is strictly not valid close to the resonance, while 
our approach is only valid around resonance because we 
ignore the negative energy contribution to the scattering 
amplitude, i.e. the terms that involve creation of vir- 
tual electron-positron pairs in the intermediate states in 
the second order perturbation theory (see Sec. II C of 
Ref . IToT) . Another difference from Blume's theory is the 
fact that our work is based on fully relativistic quantum 
mechanics, while Eq. I|38|) exploits the semi-relativistic 
approximation. This difference makes direct comparison 
of the two theories difficult. This has been discussed by 
Strange ^l| who has rederived Eq. 1)38(1 as the nonrela- 
tivistic limit of a fully relativistic theory of x-ray scat- 
tering. For these reasons and the fact that there is no 
one-to-one correspondence between the terms in our ex- 
pression for the scattering amplitude and Blume's expres- 
sion, there is no straightforward way to compare the two 
theories. It is often stated that Blume's expression Eq. 
(|38[) shows that the cross section for magnetic scatter- 
ing will yield the orbital and spin moment of a material 
separately. Although this will usually be the case it is 
not rigorously true. Eq. I|38|l cannot be applied immedi- 
ately because the initial and final states \a > and \b > 
are general many-body states that have not been defined 
in detail. For implementation purposes they must be 
described as many-electron states that will contain the 



index j which is being summed over in Eq. (|38|) . In a 
magnetic material the radial part of the basis functions 
of the single particle wavefunctions, as well as the angular 
part, depend on mi. So, we would expect the total scat- 
tering amplitude to have a contribution from the orbital 
angular momentum associated with each single-particle 
state, but this is not the same as being proportional to 
the total orbital angular momentum. For example a two- 
particle state composed of two single- particle states with 
mi = ±1 has the same z-component of orbital angu- 
lar momentum as a two-particle state composed of two 
single-particle states with mi — ±3, but Eq. (|38|l docs 
not suggest that they will have the same scattering am- 
plitude. Nonetheless, Blume's expression implies that a 
strong dependence of the cross section on the components 
of the magnetic moment is likely and indeed, this is ex- 
actly what we have found, an approximate, but by no 
means rigorous proportionality between orbital moment 
and magnitude of the cross section which is dependent 
on the polarization of the x ray. Furthermore, Figure 
demonstrates explicitly the dependence of the cross sec- 
tion on the magnitude of mi of the occupied individual 
electron states. 

The question that now arises is how our computed x- 
ray scattering results can be interpreted in terms of the 
detailed electronic structure of praseodymium. In order 
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FIG. 5: The scattering cross section at the Mrv and My edges for praseodymium. Each figure is for a different SIC configuration 
but which produces roughly the same orbital and spin moment. The thick line is for negative helicity incident x rays and the 
thin line is for positive helicity incident x rays: (a) Localising the mi = —1 and mi = +1 electrons with spin up, yields an 
orbital moment of —0.07 /j,b and a spin moment of 2.46 ^b; (b) Localising the m; = —2 and mi = +2 electrons with spin up, 
yields an orbital moment of —0.12 /ib and a spin moment of 2.47 fiB', (c) Localising the mj = —3 and mi = +3 electrons with 
spin up, yields an orbital moment of —0.05 and a spin moment of 2.52 /ib- 



to understand this we analyze the electronic structure of 
fee Pr for the cases where the orbital moment is equal to 
-4.79 hb, -0.05 fiB, and 4.69 [ib in detail. We expect the 
scattering cross section to reflect the Pr /-electron den- 
sity of states. Although the shape of the cross section 
is partially determined by the density of states (DOS) 
the total DOS changes very little when pairs of electrons 
with differing orbital moments are localized. Therefore, 
a simple interpretation of the changes in the cross sec- 
tion with the orbital moment in terms of the total DOS 
cannot be made. In relativistic theories of magnetism 
different values of total angular momentum j with the 
same z-component rrij are coupled and further decompo- 
sition has little meaning [2^ . To facilitate understanding 
of the differences in the spectra as orbital moment varies 
we show a selection of density of states curves, decom- 
posed by the azimuthal quantum number mj in Figs. Eh 
03 There are several points that should be noted about 
these pictures. 

The rrij — ±7/2 (these are pure j = 1 + 1/2 states) 
figures describe / electron states with a well-defined j 



value, while all the others show / states with two dif- 
ferent values of j (j = I + 1/2 and j = I — 1/2). In all 
the pictures except rrij = ±7/2 there are two main peaks, 
however these two peaks do not necessarily have the same 
weight. The separation of the peaks represents the spin 
and spin-orbit splitting of the individual values of rrij. 
The splitting between the unoccupied /-states is around 
0.1 Ry while the splitting between the occupied and un- 
occupied states is about 0.7 Ry. The smaller narrow 
peaks in some of these figures represent the hybridiza- 
tion of different /-states between themselves. Some of 
these densities of states are markedly broader than oth- 
ers and this is a reflection of the degree of hybridization 
with the conduction s — d electrons. 

In Fig. |S|we have chosen to apply the self-interaction 
corrections to the /-electrons which correspond to (mi = 
— 3,to s = +■=) and (mi — — 2,to s = ±|) (configuration 1 
in Table I) in the nonrelativistic limit, and this is reflected 
in the density of states having a very large and narrow 
peak at around -0.7 Rydbergs for rrij = — 5/2 and —3/2. 
There is nothing for these states to hybridize with so 
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FIG. 6: The I = 3 contribution to the density of states of praseodymium decomposed by the rrij quantum number for the 
case when the m ; = —3 and itij = — 2 / states with spin up are occupied (localized). In the top right of each figure is the 
self-consistently calculated orbital moment. Each figure is also labelled with the relevant value of the rrij quantum number. 
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FIG. 7: The I = 3 contribution to the density of states of praseodymium decomposed by the rrij quantum number for the 
case when the m ; = —3 and itij = +3 / states with spin up are occupied (localized). In the top right of each figure is the 
self-consistently calculated orbital moment. Each figure is also labelled with the relevant value of the rrij quantum number. 
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they are very tall and narrow atomic-like states. For 
nij = — 5/2 and —3/2 the density of states has a more 
band-like component corresponding to a single electronic 
state just above the Fermi energy. For most of the other 
values of rrij there is a density of states corresponding 
to two electron states close to €f and for rrij = 7/2 the 
density of states close to £f corresponds to a single pure 
j = 1 + 1/2 state. 

In Fig. U| we have selected the /-electrons which corre- 
spond to (to/ = — 3,to s = +i) and (to/ = +3,to s = +|) 
in the nonrelativistic limit for the SIC (configuration 6 in 
Table I). Here it is the rrij = — 5/2 and the rrij = +7/2 
components of the density of states that have the local- 
ized state around -0.7 Rydbergs below the Fermi energy. 
This means there is no rrij = 7/2 character around ep 
at all in this case. For most other values of rrij we can 
clearly see that there are two f-states close to e-p- De- 
tailed examination of these peaks shows that the domi- 
nant cause of the splitting is the exchange held, although 
the splitting is also influenced by the spin-orbit interac- 
tion. For rrij = — 7/2 there is only one state close to ep 
of course. In Fig. [S] we have chosen to apply the self- 
interaction corrections to the /-electrons which corre- 
spond to (to/ = +3,to s = +i) and (to/ = +2,m s = +A) 
in the nonrelativistic limit (configuration 11 in Table I). 
This time it is the rrij = 5/2 and rrij = 7/2 states that 
are localized, and again there is no rrij = 7/2 character 
around €f. The rrij = —5/2 and rrij = —3/2 state have 
the spin-split behaviour close to cf in this case. The 
other values of rrij behave as before. 

It is clear from figures [£] to |H] that in some rrij chan- 
nels there is a small amount of band-like /-character be- 
low the Fermi energy. This indicates that there are two 
types of /-electron in our calculation, the localised /- 
electrons which determine the valence and the delocalised 
/-electrons which determine the valence transitions. 0] 
It is the delocalised /-electrons that are principally re- 
sponsible for the non-integer values of the orbital mo- 
ments shown in Figures [3] and 0] (although there is also 
a small contribution from the valence s — d electrons). 

Comparison of the corresponding diagrams in Figures 
El [3 and [8] shows dramatic differences. Even though the 
total density of states is fairly insensitive to which /- 
electron states are occupied, the TO^-decomposed density 
of states is obviously drastically altered depending on 
which electrons are localized. In particular the /-states 
just above the Fermi energy form a significant number 
of the intermediate states in the formal theory described 
earlier. Therefore if key ones are localized they become 
unavailable as intermediate states for the spectroscopy 
and the cross section may be substantially altered. Of 
course, occupying one /-state means that some other / 
state is not occupied which may then also play a role 
as an intermediate state for the spectroscopy. Indeed, 
how much the unavailability of particular rrij substates 
affects the spectra depends on other factors too, includ- 
ing the El selection rules which are composed of angular 
matrix elements. Each angular matrix element contains 



four terms in the form of a product of Clebsch-Gordan 
coefficients and a geometry and polarization dependent 
factor. A further influence is the fact that the LMTO 
coefficients AL A (defined in Eq. (17) and completely de- 
termined by a self-consistent band structure calculation) 
associated with the /-electrons are found to be fairly in- 
dependent of the rare earth element under consideration 
but their magnitude has a clear but complex linear pro- 
portionality to TO/. 

Detailed analysis of the major contributions to 
the cross section suggests that the highest peak is 

formed by the core-to-valence transitions (^d^rrij^J — > 
(fs,mj + (-)l) for the M w LCP(RCP) edge scat- 
tering and ^ds,TOj^ — * (fr , rrij ; + (— )1^ for the My 

LCP(RCP) edge scattering. The former transition for 
Mrv case is in agreement with the nonrelativistic selec- 
tion rule which forbids a Aj = 2 transition, although 
this transition is not totally forbidden in the relativistic 
El selection rule. In the My case, the Aj = transi- 
tion is observed to form part of the shoulder rather than 
contributing to the main peak. Furthermore, within the 
transitions forming the main peak, the contribution to 
the LCP scattering at both the Mrv and My edge is 
the largest from the most positive allowed rrij value of 
the core state. On the other hand, the most negative 
rrij value of the core state gives the largest contribution 
to the RCP scattering. This indicates the fact that the 
Clebsch-Gordan coefficients which are used to calculate 
the selection rules are a dominant factor in determining 
the relative size of the cross section peaks. The origin 
of this is simply in the properties of the Clebsch-Gordan 
coefficients which vary smoothly between either and 1 
or and -1 depending on the values of the other quantum 
numbers. 

From these considerations, we see that the separation 
of the LCP and RCP peaks by 1 to 2 electronvolts is a 
reflection of the spin-splitting of the states. In relativistic 
theory m s and to/ are not good quantum numbers. Fur- 
thermore, because of the magnetism, different values of 
j with the same rrij are also coupled. However, it is still 
possible to associate (s z ), (l z ) with these quantum num- 
bers and also to recognize the dominant j in atomic-like 
unhybridized bands. For example, in the case of LCP 
scattering at the Miv edge, the largest contribution to 
the cross section comes from (k = 3, rrij = + §)-like or- 
bitals. The two 4/ states which have this rrij as the main 
contributor are characterized by ((s z ) — +h, (l z ) — +2) 
and ((s z ) ~ — A, (l z ) ~ +3). Electronic structure cal- 
culation shows that the former state is dominated by 
k = —4 and the latter by k = 3. Therefore the Miv 
LCP peak is most affected by the availability of the spin- 
down (l z ) ~ 3 state as an intermediate state. Similar 
analysis shows that the Miv RCP peak is most affected 
by spin-up (l z ) ~ —3 state, My LCP by spin-up (l z ) ~ 3, 
and My RCP by spin-down (l z ) ~ —3 state. 
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FIG. 8: The I = 3 contribution to the density of states of praseodymium decomposed by the rrij quantum number for the 
case when the m ; = +3 and itij = +2 / states with spin up are occupied (localized). In the top right of each figure is the 
self-consistently calculated orbital moment. Each figure is also labelled with the relevant value of the rrij quantum number. 
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Although this analysis is a gross simplification, it does 
explain why the relative peak energy positions in the LCP 
and RCP scattering cases swap between the Mrv and My 
edges (see Fig. |SJ. Of course this is true only if these 
states are still available after the chosen localizations by 
SIC. The effect of localization on the MXRS spectrum 
is most dramatic if SIC is applied to these key states, 
changing the peak energy separation as well as the scat- 
tering amplitude between the LCP and RCP scattering 
cases. 

Some empty valence band / states participating in the 
scattering process have nearly equal mixture of the two 
j characters, i.e. j = 1 + 1/2 and j = I — 1/2. If there 
is strong spin-up and spin-down character in the unoccu- 
pied valence states described by a specific rrij then both 
spin states may be available as the intermediate states 
for the spectroscopy. Thus we may clearly see a two-peak 
structure in the rrij decomposed amplitude for a certain 
polarization at the absorption edge. Figure El shows the 
core rrij decomposed LCP scattering amplitude and the 
two-peak structure mentioned above is clearly visible for 
rrij = +i at the M IV edge. 

In certain cases we can interpret the apparent re- 
lation between the magnetic cross section and the z- 
component of the total orbital moment as follows. Be- 
cause (l z ) + (s z ) = (j z ) holds, then we see that if we 
apply self interaction corrections to states systematically 
according to Hund's rules, what is effectively done is to 
occupy the states in order of rrij. As stated earlier the 
rrij decomposed relativistic magnetic scattering cross sec- 
tion has a 'proportionality' to rrij due to the Clebsch- 
Gordan coefficient in the angular matrix element expres- 
sion defining the El selection rules. Whether this pro- 
portionality is direct or inverse depends on the polar- 
ization of x rays. In addition, according to the electronic 
structure calculation, as the unhybridized state goes from 
((l z ) ~ -3, (s z ) ~ |) to ((I,) ~ +3, (s z ) ~ i), the domi- 
nant j changes from j = | to j = | gradually. This tells 
us two things. Firstly we notice that if a certain state has 
a major impact on the scattering cross section at the -Mrv 
edge for RCP photons, then this same state has a rela- 
tively minor effect on the cross section for LCP photons 
at the same edge because of the Clebsch-Gordan factor 
in the expression for the El selection rules as mentioned 
above. Secondly we see that this same state also has 
only a minor effect on the My cross section because the 
value of j for the intermediate states involved in major 
transition differ between Mrv and My. 

As the SIC configuration varies from ((Z z ) ~ —3, (s z ) ~ 
+ i) and ((/,) ~ -2,{s z ) ~ +i) to {{l z ) ~ -3, (a.) ~ 
+ \) and ((l z ) ~ +3, (s z ) ~ +^) so that there is a sys- 
tematic change in the z-component of the total orbital 
moment, the Mrv RCP cross section increases because 
the second, third and so on, strongest contributors to the 
cross section become additionally available as intermedi- 
ate states as they are released from the SIC localization. 
However, they have progressively less impact as we pro- 
ceed through this series of quantum numbers since the 



major j gradually changes to j = |. The cross section at 
the Miv edge for LCP photons is not affected much by 
this change in quantum numbers since neither the initial 
nor the final SIC combination in the above series involves 
the major contributors to Mrv LCP cross section. On the 
other hand, the My edge LCP cross section is reduced as 
more and more significant contributors are removed from 
the available intermediate states, while the cross section 
at the My edge for RCP photons is not much affected for 
the same reason as Myy LCP case. Obviously the above 
change in SIC configuration is very artificial. However 
as the states are filled up according to Hund's rule as 
we proceed through the rare earth series, we would ex- 
pect to observe changes in the cross section governed by 
these considerations for rare earths where the intermedi- 
ate states can be considered as atomic-like. However, a 
very different interpretation of the x-ray spectra may be 
required in the case where dclocalizcd band-like interme- 
diate states are of primary importance, as is the case in 
resonant x-ray scattering at the K and Lnjn edges. 

Finally, we are unaware of any experimental measure- 
ments of the MXRS spectra of praseodymium or it com- 
pounds at the Mjy or My edge. However a careful com- 
bined neutron [53 and x-ray (at the £11,111 edges) 
investigation into the magnetism of Ho^Pr^^ alloys has 
concluded that the Pr ion does have a 4/ moment at 
all values of x. Deen et al. |5{| have performed MXRS 
measurements at the L edges in Nd/Pr superlattices and 
found a large peak at the absorption edge and a high en- 
ergy shoulder corresponding to dipolar transitions to the 
broad 5d band. We hope that our calculations will stim- 
ulate detailed experimental x-ray studies of Mrv and My 
edges of Pr, in pure Pr and in its alloys and compounds. 

V. CONCLUSIONS 

In conclusion, a theory of magnetic x-ray scattering 
that is based on the LSD with self-interaction correc- 
tions and second order time-dependent perturbation the- 
ory has been described. We have illustrated the theory 
with an application to fee praseodymium and used this 
example to illustrate the dependence of the scattering 
cross section on spin and orbital magnetic moments. It 
has been shown that the theory quantitatively reproduces 
the dependence on the spin and orbital magnetic mo- 
ments originally predicted qualitatively . 
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